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Abstract 

We consider anisotropic colloidal particles immersed in a solution of long, flexible, and nonad- 
sorbing polymers. For the dumbbell shapes of recently synthesized particles consisting of two inter- 
secting spheres and for lens-shaped particles with spherical surfaces we calculate the isotropic and 
anisotropic interaction parameters that determine the immersion free energy and the orientation- 
dependent depletion interaction between particles that are induced by the polymers. Exact results 
are obtained for random-walk like (ideal) polymer chains. 
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I. INTRODUCTION 



In colloidal suspensions containing polymer chains, there are tunable effective interactions 
between the colloid particles. Free nonadsorbing polymer chains avoid the space between 
two particles, leading to an unbalanced pressure, which pushes them towards each other. 
Such depletion forces for an isolated pair of immersed particles or for a single immersed 
particle near a wall have been measured in recent experiments^. 

Here we consider anisotropic colloid particles. For the dumbbell shapes of recently 
synthesizecl^ particles consisting of two intersecting spheres, and for lens-shaped parti- 
cles with spherical surfaces, as in Fig. 1, we calculate the immersion free energy and the 
orientation- dependent depletion interaction. The predictions are compared with results for 
prolate and oblate ellipsoidal^, which have also a symmetry axis of revolution and a sym- 
metry center of reflection. 

The case of large particle to polymer size ratio can be investigated by means of small 
curvature expansions of the Helfrich or Derjaguin type, but here we consider mesoscopic 
particles which are small compared to characteristic polymer lengths such as the gyration 
radius TZg, and we concentrate on the case of ideal, random- walk like, polymers. The well 
known correspondence^ between the statistics of long flexible polymers and critical field 
theories allows us to use the small particle operator- expansiorJ'^'^ for predicting the polymer- 
induced interactions. The operator weights in the expansion for dumbbells and lenses are 
calculated by a conformal mapping to a wedge-geometry. 

We introduce the polymer-magnet analogy and small particle expansion in Sec. II, discuss 
density-profiles in a wedge and outside a dumbbell or lens in Sec III, and evaluate, for ideal 
polymers, the corresponding small particle amplitudes in the Gaussian model in Sec IV. 
These results are compared with corresponding results for ellipsoids and, in Sec. V, with 
a more general class of weakly anisotropic particles. In Sec. VI the amplitudes are used 
to determine the orientation-dependent interactions, and in Sec VII we summarize the new 
results. Some technical details are relegated to Appendices A-D. 
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II. POLYMER-MAGNET ANALOGY AND SMALL PARTICLE EXPANSION 



In the polymer-magnet analogy the partition function of a polymer chain with ends at ri 
and r2 corresponds to the order-parameter correlation function {'^12) of a Ginzburg-Landau 
model or field theorj^. Here ipi2 is the product $(ri) $(r2) of two order parameter fields 
$. Ideal polymers correspond to a Gaussian Ginzburg-Landau model with Hamiltonian 

"1 o t 



H = dr 



(2.1) 



where the integration extends over the volume outside the particles, and where the order 
parameter satisfies the Dirichlet boundary condition 

$ = (2.2) 

at the particle surfaces, since we consider nonadsorbing polymers. We always consider length 
scales much larger than the persistence and extrapolation lengths. 

The free energy F it costs to immerse particles in a dilute polymer solution in an un- 
bounded space or in the half space bounded by a wall is determined by the polymer partition 
functions with and without the particles and is given by^ 

F/po = -/^ y dridr2[{(pi2)n+sn- {'fi2)n]- (2-3) 

Here po = nksT is the ideal gas pressure in the dilute solution with chain density n, and H + 
STi and Ti denote Ginzburg-Landau Hamiltonians of the form ()2.1|) in presence and absence 
of the particles, respectively. The dependence of the double integral on the temperature 
deviation t from the critical point is converted into the dependence on TZ^ of —F/pq by 
means of the inverse Laplace transform £(..) = J {dt/ 27! i)exp{3t7l'^/d) (..) where d denotes 
the spatial dimension. 

Consider particles with a size much smaller than TZg and a shape that is symmetric about 
both a center of reflection and an axis of revolution. Examples are rods, disks, ellipsoids, 
dumbbells, and lenses. In the spirit of the operator-product expansion, a small mesoscopic 
perturbation in a critical field theory can be represented by a sum of point operators. Thus, 
for a single small particl o^i^i^i^° with center at rp. 



e"^^ oc 1 + a/ + (T^ , (2.4) 
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where 



(Ti = ai e(rp) + ... , (2.5) 

aA = &i9je(rp) + 62r||||(rp) + ... (2.6) 

are linear combinations of isotropic(I) and anisotropic(A) operators from the operator alge- 
bra of the Ginzburg-Landau model, which reflect the symmetries of the particle shape and 
boundary condition. Here e oc — $^ is the energy density, d\\ is a derivative along the particle 
axis, and T|||| is the diagonal component of the stress tensor of the field theory along the 
axis. 

Only the operators of lowest scaling dimensions x = d — l/u, d + 2 — l/u, d are shown in 
Eqs. ()2.4|) - ()2.6p . and their coefficients ai, bi, 62 by scale invariance must be proportional to 
the particle size raised to the power x. For ideal chains (Gaussian model) the Flory (correla- 
tion length) exponent u equals 1/2, so that (9ye and T|||| have the same scaling dimension d, 
and the 61 and 62 terms both contribute to the leading anisotropic behavior. For chains with 
excluded volume interaction (corresponding^ to the A^- vector model in the limit ^ 0), u 
is larger than 1/2, and we expect that the 62 term dominates the anisotropic behavior of a 
small particle. 

The coefficients ai, 61, 62 depend on the size and shape of the particle but are independent 
of other distant particles, of the distant boundary wall of the half space, and^^ of t . Thus 
we evaluate the coefficients for a single particle in an unbounded space at t = 0, and then 
use them to make predictions for the interaction between particles or a particle and a wall. 
The coefficients can be evaluated from the density profiles (e(r)) and {Tki{r)) of the energy 
density and the stress tensor that are induced by a single particle. While ellipsoids have 
been considered in Refs. Oj^ we concentrate here on dumbbells and lens-shaped particles. 

III. DENSITIES IN A WEDGE AND OUTSIDE A DUMBBELL OR LENS 

A system at the critical point containing a colloidal dumbbell composed of two overlapping 
spheres, or a lens with two spherical surfaces, can be conformally mapped onto a critical 
system filling a wedge with opening angle a, which is smaller or larger than vr in case of the 
dumbbell or lens, as in Fig. 1. As explained in more detail in Ref. 0, an inversion about 
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the point denoted by the heavy dot on the left hand side of Fig. 1 maps the interior of 
the wedge onto the exterior of a particle with a dumbbell or lens shape. The two boundary 
half planes of the wedge and the edge where they meet are mapped onto the two spherical 
surfaces of the particle and the circle C of diameter V where they intersect. V is related 
to the diameter L of the two spheres by P = Lsin(a/2). For a = 0, vr, and 2-7? the particle 
on the right hand side of Fig. 1 becomes a dumbbell of two touching spheres, a spherical 
particle, and a circular disk, respectively. 

For the wedge the boundary-induced density profile of a scalar operator O, such as the 
energy density e, has the fornJ^J^iii 

(0(fe,p,l^))tl = B]i'p~-op^,{a,n) . (3.1) 

Here the position vector r is expressed in cylindrical coordinates (re,p, fi), where the edge 
of the wedge is the axis. The component is parallel to the edge (and in general has 
dimension d — 2), and the two component vector perpendicular to the edge is determined by 
its angle Q with the symmetry half plane of the wedge, i.e. —a/2 < ^2 < «/2, and by the 
distance p from the edge. The indices characterize the surface universality classesi^ of 
the two boundary half planes Q = {—a/2, a/2) of the wedge. While P is a universal scaling 
function, Bq is the non-universal amplitude in the pair correlation function 

(O(r)0(O))buik = 5or-2-« (3.2) 

of O in the bulk. 

The corresponding density (0(r))particie outside a dumbbell or lens follows^^ from 
and the conformal mapping and has the form given in Eqs. (5.14) and (5.10) of Ref. 
Turning to the energy density O = e with scaling function V = B and equal boundaries 
i = i, the behavior 

4,(«, VL) = et\a)[l + ef{a) Vi^ /2 + ...] (3.3) 

near the symmetry half plane f2 = of the wedge determines the profile (e(r)) particle far from 
the dumbbell or lens and yields^ 

af =V''^ef{a)/B\l'' (3.4) 



and 



for two leading coefficients in the small particle expansion. 

The boundary-induced density profile (TKA(r)) wedge of the stress tensor in the wedge is 
given by 

{T.x{h,P,n)tg. = p-'nA^)[S.x-du^:\n)uPm, (3.6) 

where u*^"'^(ri) is the unit vector normal to the half plane f2 = const which contains r, as 
in the left hand side of Fig. 2. Unlike p.ip there is no nonuniversal amplitude in ()3.6p . 
and the p-exponent and the f2-dependence are trivial. Only the variation of the universal 
amplitude r with the opening angle a of the wedge depends on the bulk universality class 
and the surface classes i,j, and remains to be determined. The stress tensor density ()3.6p 
has a vanishing trace and obeys the continuity equation, as discussed in Appendix A. 

We note two special cases: (i) For a ^ it and i = j the wedge becomes the half space 
with a uniform boundary, and r —>■ since the stress tensor density vanishes^^. According 
to the left hand side of Fig. 2, only for r = is ()3.6|) consistent with the symmetries of the 
half space, (ii) For a — >• 0, r diverges as 

r,, a-" (-A,,) , (3.7) 

where Aj j is the universal amplitude which determines the stress tensor density in the 
parallel plate geometryi^. If the width of the ffim is u and both tensor components are 

parallel to the plates, (Tparallel,parallel)film = 

Using the inversion transformation for the conformal stress tensoi'2ii& one finds from Eq. 
(I3.(i|l the stress tensor density 

{Tkiir))^;^Le = (V/AYnAo^) [hi - dui''\r)u\''\r)] (3.8) 

outside the dumbbell or lens. Here 

A' = ^[r2 - (p/2)2]2 + I)V2 (3.9) 

with r and ry the distance of point r from the particle center and its component parallel 
to the particle rotation axis. The inverse length T) / in ()3.8|1 equals 6(r)/p(r), where h 
is the dilatation factor |det(9r/9r)|^/'^ of the conformal mapping, and diverges as {r\\^r) — >■ 
(0, P/2), as the circle C of intersection is approached. The unit vector u(^) in the particle 
geometry is the counterpart of u^"^ in the wedge geometry. It points along the surface normal 
at r of the spherical surface portion Sc,r which contains r and is bounded by the circle C. 
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Sc,r is the image of the half plane = const which contains f . With the particle 
one of the Cartesian directions, 



u 



r) = {4||[r2 - {V/2f] - 2r,^} . (3.10) 



The vector field u*^^^ for given T) is shown on the right hand side of Fig. 2. Both u*^^) and 
u'-"'' are independent of a and the bulk and surface universality classes 
In leading order P -C r, 

{TM{r))^;;^Le - (V/rYr^Ac^) [Sm - dl,,nlm] , (3.11) 

with 

1k,m = 1k,mij^) = 5k,m — '^fk^'m/ ^'^ ■ (3.12) 

On comparing with the stress tensor correlation function 

(rfc,(r)r^„(0))bulk = BTr'^\{l/2)[Xk^^ll^r^+Tu,n1l,m] ' iM d)5ui5mn} (3-13) 

in unbounded bulk^^^^^, Eq. fl3.1ip implies that 

(7^;t/(r))!i;lie - feF(r/cKr)7||||(0))buik , l^«r, (3.14) 
with the stress tensor contribution in the small particle expansion given by 

6j'^')T|||| = -P%(a)^T||||. (3.15) 

JDnp 

For the special case a ^ 0, where V aL/2 and Eq. ()3.7j) applies, Eq. ()3.15|1 reduces to 
the expression 63 '•'^ {L/2Y^jd/ Bt for a dumbbell of two touching spheres with diameter 
L, given in Eq. (2.15) of Ref. 
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18 



The form of (2l)wedge f|3-(jp has been calculated for special cases. See Refs. 
for d = 2, Ref. llJ for symmetry breaking surfaces i = j as d ^ A, and Appendix B for 
Dirichlet boundaries and a = tt/2. Here we indicate how ()3.(i|l can be derived in the general 
case. For r on the particle rotation axis, {T)^^2^^Ye niust have an eigenvector parallel to 
the axis, by rotation symmetry. That the particle axis is the image of a circle in the plane 
f e = of the wedge with center in the edge and passing through the center of inversion, 
see the long dashes in Figs. 1 and 2, implies an eigenvector u*^") of {T)^'^J^^ tangent to this 
circle. Likewise, degenerate eigenvectors perpendicular to the axis imply d — 1 degenerate 
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eigenvectors perpendicular to u"^"). The simple form (j3.6j) of the symmetric, traceless, and 
conformal tensor density (T) for the wedge then follows from scaling (dilatation symmetry), 
translation, and reflection symmetry in the edge-subspace, and the f2-independence of r is 
due to the continuity equation of the stress tensor, see Appendix A. In Appendix B we also 
discuss, within the Gaussian model, the more complicated form of the density (r*^"^^'^))wedge 
of the canonical stress tensor, which also obeys the continuity equation, but is not trace-free 
and not a conformal tensor, and for which u(")(i7), in general, is not an eigenvector. Its 
eigenvalues depend on both p and Q, and some of the eigenvalues diverge as the boundary 
planes of the wedge are approached. The simple form ()3.6|) is recovered on adding the 
'improvement-term—. 

We briefly comment on particle shapes where the horizontal axis on the right hand side of 
Fig. 1 (passing through the small triangle and square) is the rotation axis. These particles 
resemble an apple (self-intersecting torus) for a < tt and an american football for a > vr and 
could be conformally mapped onto a cone with opening angle a. 

IV. GAUSSIAN MODEL WITH DIRICHLET BOUNDARIES 

In d = 2 spatial dimensions, a wedge can be obtained from the half plane by means of 
an appropriate conformal transformation. The dependence on the angle a only enters via 
the transformation and is to a large extent model-independent. For example, for i = j, the 
stress-amplitude r is independent of the surface universality class i and equals [(vr/a)^ — 
l]c/(247r), where only the universal bulk-constant c, the 'conformal charge', depends on the 
model class. Using ()3.7|) . this is consistent with the i-independent film-amplitude Aj j = 
— 7rc/24 for a strip with equal boundary conditions'^. Also the form of the scaling functions 
Vi^i = A^(!j)[{a/'K)cos{7T^l/a)]~^° is to a large extent model-independenlii^ and is completely 
determined by the bulk scaling index xq and the universal, but z-dependent, half-space 
boundary amplitude A^^ of the scalar operator O. For i not equal to j, the stress amplitude 
Tij = — (c/(247r)) — Ajj/a^ also depends on the boundary universality classes. Since Ajj 
for i 7^ j is in general different from — 7rc/24, Tij is nonvanishing even for a = tt, i.e. for the 
half plane with a non-homogeneous boundary^. 

In (i > 2, however, no such conformal transformation exists, and the a-dependent quanti- 
ties V and T depend on the bulk and surface universality classes in a much stronger way. In 
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particular a bulk amplitude (like c) and the parallel plate amplitude Ajj, are not sufficient 
to predict Tij{a). 

For the Gaussian Ginzburg-Landau field theory (|2.ip at the critical point t = in d > 2 
spatial dimensions, the scaling function S of the energy density in the wedge is given by 

^ = -($2(f))„edgep'-V(V2^d) , ^rf = 7r-^/'r((rf/2)-l)/4, (4.1) 

where we have suppressed the indices = (D, D) with D for Dirichlet^. The stress 

tensor in the Gaussian model is the sunt^^ 

T.x{^) = TirH?) - hx{^) (4.2) 

of the canonical stress tensor 

Tir\r) = {d.m)0xm) - s.,{vmr/2 (4.3) 

and the so-called improvement term with 

4A(i-) = - 5.AA]$^(f) . (4.4) 

Unlike T^^^^\ the density of T in (j4.2j) has a vanishing trace and transforms as a conformal 
tensor, see Appendices B and C. The bulk amplitude in (|3.13p equals [Sd{d — 2)]'^d/{d—l) 
for the Gaussian model. 

The boundary-induced densities ($^(r))„cdgc and (T«;A(r)) wedge follow in an obvious way 
from the boundary-induced contribution 5(<^i2)wedge of the propagator in the wedge 

(<^12)wedge = ('^12)bulk + '^(<^12)wedge , (4.5) 

with 

0,2 = HhMh) . (4.6) 

For the special opening angles a = n/g, with g a positive integer, 5((^i2)wedge is a linear 
combination of 2(7 — 1 bulk propagators from ri to images of r2, see Appendix B. 

In order to determine Cq, 62 in (|3.3|) and r in (|3.6|) for arbitrary a, we use the representation 



of Ref. 



20I at the critical point 

» 00 
(^i2)wedge = J d''~'q{2nf-' e^'i^^--^-) (2/«) Im./MP<)K„../MP>) 

m=l 

r /I Oo \ 1 

(4.7) 
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a /. 



where p< = min(pi,p2), P> = niax(pi,p2), and / and K are modified Bessel functions. A 
more explicit expression arises for (rei,pi) = (re2,P2)- Using 



r dxx'^-'l^^/MKm./a{x) = 2^-^- %./ \. f dtt^-'^t^-'{l-t)-^^ (4.8) 

Jo 1 (1 — A) JO 

for 2A = d — 2 in order to rewrite the g-integral in a form where the m-summation can be 
done^, one finds 

^d^^i2)wedge|(f,i,pi)=(fe2,P2) = -(2/a)sin(7rrf/2)p2"'^ / rft^d(t) 

J 

J 1 - W"cos[(r]i - ^2)71 /a] 1 + W°cos[(l]i + ^2)7^/ a] \ 

^ \ 1 + - 2tV"cos[(l]i - fi2)7r/«] 1 + tW" + 2W"cos[(l]i + ^2)7x1 a] J " 

where p = pi = p2 and 

= t^'^-^^/^ ^1 - (4.10) 

The expression in Eq. (14. 9|) satisfies the Dirichlet boundary conditions since the curly bracket 
vanishes for Qi (or ^^2) equal to a/2 or —a/2. For Qi ^ Q2 the t-integral converges in the 
interval 2 — (27r/a;) < d < 4 of spatial dimensions d with the lower and upper limits coming 
from the behaviors of the integrand for t \ and t y 1, respectively. The bulk divergence 
{p\fli — fl2\Y~'^ for \Qi — ^2! and d > 2 contained in ()4.9j) comes from integrating the 
first term in curly brackets over a region near t = 1 where 1 — t is of the order of — i72|- 
A convenient way to evaluate the boundary-induced density 

(<l>^(r))„edge = limr^_,f,f2_,f (5((^i2)wedge (4.11) 

and its scaling function S{a,Q) in 1)4.11) for = 3 by means of Eq. ()4.9|) . is based on the 
observational that S is an analytic function of d in an interval that includes both d = 2 and 
(i = 3 as interior points. Since the limit ri r, r2 — > r of the bulk propagator oc |ri — r2p^'^ 
vanishes for d < 2 (while it is infinite for d > 2), one may replace 5(<^i2)wedgc by ((^i2)wedge 
and use ()4.9p . in calculating the limit in ()4.1H) for d < 2. This leads to 



£{a,Q) = {\^/a)sm{7id/2) f\t^d{t) 

Jo 

j 1 1 + W"cos(27rfi/a) 



d<2 , (4.12) 



[1-W" l + t27r/" + 2W"cos(27rfi/a) j 

in terms of an integral which is well defined for 2 — (2,71 /a) < d < 2 and which has to be 
analytically continued^i in order to obtain S for d = 3. 
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For Co (a) = £{oi, Vt = Q) the integral in ()4.12|) becomes 

J(a, d) = C dt m,{t) -4^ , (4.13) 



(4.14) 



1 _ t'^-K/a 

and the continuation can be made by rewriting J as the sum of J^^^ and J*^^-*, where 

2 



j(i)(a,rf) = rft^d(t)r/" 

JO 



-lit) 



and 



j(2)(ct^rf) = dt^d{t)t'''''l{t). (4.15) 

JO 



Here 



= ^ + (4.16) 

are the first two terms in the Laurent series of 2/(1 — t^'^l"') around t = 1 so that the 
integrabihty domain 2 — (27r/«) < < 4 of J'-^-' extends up to c? = 4. As a sum of beta- 
functions the continuation to d = 3 of the integral J*^^-* is trivial and yields J^'^\a^ 3) = — 1, 
implying 

eo(a) = (V2/a)[l - ^^^^",3)] , ci = 3. (4.17) 

Here J*^^-'(q;,3) follows from the right hand side in Eq. ()4.14|) on replacing by \l/3 = 
t^^l'^ I (1 — t), and for a arbitrary between and 27r we have calculated it numerically. Using 
Eq. ()3.4p . the corresponding results for aiB]^"^ /L = sin(a/2)eo(«) in the interval < a < vr 
and aiBj/yV = eo(a) in the interval vr < a < 27r are shown in Fig. 3. Analytical results 
for some special values of a are given in Table 1. 

In the Gaussian model some amplitudes of isotropic operators beyond leading order in 
()2.5|) . such as $^ and are also determined^S by ai. 

To calculate the small particle anisotropy-amplitude bi in ()3.5|1 . we need the coefficient 
oc of S. The contribution of order Q"^ to the curly bracket in ()4.12j) contains a factor 
1 — f^^" and leads to a convergent integral up to d = 4. Thus, no continuation is necessary, 
and 

eo(a)e2(a) = (4y27rVa^) dt*3(t)r/"(l-r/")/(l + r/")^ , d = 3, (4.18) 

JO 

which on using ()3.5p with ()3.4|) . ()4.17|) leads to the results for bi shown in Fig. 4 and in 
Table 1. 
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In Appendix D we use similar continuations in d to calculate the stress tensor amplitude 
r(a) from ()3.6|) and ()4.2|) - ()4.4|) in ci = 3 with the result ()D15|) . The second anisotropy- 
amplitude 62 then follows from (|3.15p and the value d/Bx = 327r^, with the results shown 
in Fig. 5 and Table 1. 

It is interesting to compare dumbbells and lenses with ellipsoids. We compare a dumbbell 
with a between and vr with a prolate ellipsoid that circumscribes the dumbbell, touches it 
at the highest and lowest points, and has the same curvature at these points. Denoting by 
D\\ and D± the diameters of the ellipsoid parallel and perpendicular to the rotation axis, 

i^ll = 2Lcos2(a/4) , = V2Lcos{a/A) . (4.19) 

Similarly, we compare a lens with a between vr and 2ti with an oblate circumscribing ellipsoid, 
so that 

D\\ = Vctg{a/4:) , = V (4.20) 

where X'ctg(Q;/4) is the width^ of the lens. The amplitudes ai, 61, 62 of the ellipsoids are 



ll||2a and Eqs. (180)-(185) in Ref. 



y 



shown as circles in Figs. 3-5. They follow from Refs 
where the long axis [D\\ , and short axis , D\\\ of the [prolate , oblate] ellipsoid is 
denoted by I and s, respectively. 

As expected, the isotropic and anisotropic perturbations of the polymer system from 
dumbbells are weaker and stronger, respectively, than from the circumscribing prolate el- 
lipsoids. The oblate eUipsoids have stronger isotropic-perturbation amplitudes ai and also 
slightly stronger anisotropic amplitudes hi and 62 than the lens. 



V. WEAK ANISOTROPY 



Consider the amplitudes oi, 61, and 62 for particles with a surface S' which deviates only 
slightly from the surface 5* of a sphere with radius R. S' is obtained by shifting each surface 
point vs of 5 by a small amount ri{9s) toward the center of S at the origin. Here 9s is 
the angle which encloses with the particle rotation axis, and we consider particles with a 
center of reflection so that rj{6s) = ri{n — 9s)- For a dumbbell or lens with a = tc + 6a, we 
choose R = L/2 and obtain 

r] = 5a {L/ A) \cos9s\ (5.1) 
12 



to first order in 5a. As expected from Fig. 1, 77 in ()5.ip is non- analytic at 6*5 = 7r/2. At tlie 
end of Sec. Ill we liave introduced apple- and football-sliaped particles with a = t[ + 5a 
smaller and larger than tt, respectively. Nearly spherical particles of this family are generated 
by 

7] = 5a (L/4) siiiOs ■ (5.2) 

For weakly anisotropic prolate and oblate ellipsoids with D\\ > D± and D\\ < D±, respec- 
tively, 

D\\-D^ . , , 

V = 2 ^^"^^^s, (5.3) 

if we choose R = D\\/2. 

In the presence of the weakly anisotropic particle the propagator is given bj*^ 

{(^12) = (v?12)sphore + J dS T]{es) {T_L ±{r s) (fiu) sphere (5.4) 

to first order in rj. Here / dS is an integral over the surface 5* of the sphere, and T±± is the 
diagonal component of the stress tensor perpendicular to S. Due to the Dirichlet boundary 
condition and the form ()4.2|) - ()4.4|) of the stress tensor, Tj_^(r5') can be replaced by (9j_$)^/2 
in the correlation function in ()5.4|) with the result 

and, from the behavior of (|5.4p for R <^ ^i, r2, one finds^ 

,rf-2/,/7T (rf-2)r(rf/2) r ,,-2, 



V27rr(((i - l)/2) Jo 



and 

b2 = R""' r{{d+l)/2) Jo "^^^ ^^i^^^) [dicosOsf - 1] viOs) . (5.8) 

Explicit expressions for dumbbells or lenses, apples or footballs, and ellipsoids in an arbitrary 
spatial dimension d follow on inserting t] from ()5.1|) - ()5.3p . In = 3, 

.,y5:/L = ^-fo(^.^), 



TC IT 



62/L^ = 5a(-,--) (5.9) 
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for (dumbbell or lens, apple or football), and 

ai\ 



2v^ 3v^ ' 

&2 = {D^-D\{)D^'^-^ (5.10) 
15 

for ellipsoids. In our first order calculation D may be either D\\ or Dj^. Note that 62 is 
negative for the prolate shapes (dumbbell with 5a < 0, football with 5a > 0, and prolate 

ellipsoid with D\\ > Dj_) and positive for the oblate shapes (lens with 5a > 0, apple with 
5a < 0, and oblate ellipsoid with D\\ < D±). 

The amplitudes 61 follow from the amplitudes 62 above via the ratio 

b,^Jb, = (5.11) 
o7rv2 

For a nearly spherical particle with rotation axis and reflection symmetry, the ratio of 
the anisotropy amplitudes 61 and 62 is independent of its shape, see ()5.7j) and ()5.11|1 . but the 
ratio becomes shape-dependent for larger deviations from spherical. For example, 

I f 7 8 + V2 5 ln(2D||/Dj 



-I^.^l^j (5.12) 

for a dumbbell of two touching spheres, a dumbbell with a = tt/2, a disk, and an ellipsoidal 
needle with D\\ ^ D±. For the dumbbell-lens and ellipsoid families the modulus of the 
negative ratio biy/Bl/b2 is monotonically decreasing on increasing a from to 27r and D±/D\\ 
from to cxD, respectively, i.e. on changing from more prolate to more oblate shapes. 



VI. INDUCED INTERACTIONS 

The expressions ()2.3j) for the free energy cost F and ()2.4j) for the Boltzmann factor 
exp(— (57i) of a small particle determine the polymer- induced orient at ion- dependent interac- 
tions. In particular, the leading anisotropic interactions between a particle and a wall, 

F^lilJ = Po{cos'^p){h^/BAnV2M';{y) 

-&2[/o(y) + (l/2)/o(y/2)]}, (6.1) 

with the second derivative 



M'; = 4[/o(i/) - (l/2)/o(y/2)] 
14 



(6.2) 



of the bulk-normalized polymer density in the half spaceS- with respect to 

y = zp/TZg , zp = particle wall distance , (6.3) 
and between two particles P, Q, 

^^i's? = -Po(cosVp + cosVq) V2ai\/^/5vii?^PQ 

x{12x-^f2{x/2) + 6x-'f,{x/2) + fo{x/2)} , (6.4) 



with 



and 



X = rpQ/TZg , (6.5) 
rpQ = |rp — fqI = particle particle distance (6.6) 

/3vii = [16v^7r6i/B^ + 62] /2 , (6.7) 



follow from ()2.3|) with the half space perturbed by exp(— 57"^) oc cr^(P) and the bulk per- 
turbed by exp(— 57i) oc cr/(P)(j^(Q) + cr/(Q)cryi(P), respectively. Here /„ = i^erfc is the 
n-fold iterated complementary error function, 'dp is the angle between the particle axis and 
the surface normal of the boundary wall, and (pp,(pQ are the angles between the axes of 
particles P, Q and the distance vector rp — tq of the two particles. -F^niso^ is proportional to 
the anisotropic part of {djp + d^^Q)K{rpQ), with K the density-density correlation function 
of ideal polymers in bulk solution. 

For a dumbbell or lens with any a the particle-wall expression predicts that, for small y, 
the particle orientation parallel to the wall and, for large y, the perpendicular orientation, 
have the lowest free energies. Note that A^h has a point of inflection, and Ai'-^ is positive 
and negative for small and large y, respectively. With the values bi and 62 from Figures 4 
and 5 above, the 61 and 62 contributions both favor the same, parallel orientation for small 
y. For large y they favor different orientations, and their sum is oc — /3vii/o(2//2), with /3vii 
from ()6.7|) in which bi dominates. 

The two-particle expression predicts that particles align parallel to their distance vector, 
as expected from the attractive nature of the depletion interaction in a dilute polymer 
solution. 



Qualitatively similar 
can be taken from Refs. 



behavior applies for prolate and oblate ellipsoids, where 61 and 62 
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VII. SUMMARY AND CONCLUDING REMARKS 



We have studied the interaction between long flexible nonadsorbing polymers and meso- 
scopic colloidal dumbbells and lenses. The shape of the colloids is characterized by a pa- 
rameter a, as shown in Fig. 1 and ranges from two touching spheres for a = 0, to a sphere 
for a = TT, and to a disk for a = 27r. 

For small colloids and ideal polymers the amplitudes ai, bi, and 62 in the small particle 
expansion ()2.4|) . which determine the isotropic and anisotropic features of the interaction, 
are evaluated exactly for arbitrary a. They follow via the general relations ()3.4j) . ()3.5|) . and 
(|3.15|) from the results (j4.17j) . (j4.18p . and ()D15|) for density profiles of the Gaussian model in 
a wedge with opening angle a and Dirichlet boundary conditions and are shown as crosses in 
Figs. 3-5. We compare with corresponding amplitudes for ellipsoids that circumscribe and 
touch the dumbbells and lenses, see Eqs. fj4.19|) . ()4.20|) . Their values are shown as circles in 
Figs. 3-5. Analytical results for some special values of a are collected in Table I. 

We also consider weakly anisotropic particles of general shapes with rotation axis and 
reflection center, see Eqs. (j5.6|) - (j5.8|) . We find that the ratio 61/62 in Eqs. (j5.7|) and (j5.11|) of 
the two anisotropy-amplitudes of these particles is independent of their shape. In particular 
we consider in Eqs. fj5.9|) the shapes of a self-intersecting torus which resembles an apple 
and of an american football. 

How to obtain from the amplitudes Oi, bi, and 62 the orientation-dependent polymer- 
induced interaction between particles is discussed in Sec. VI. While the preferential align- 
ment of two identical small particles is always parallel to their distance vector, see Eq. (j6.4j) . 
the alignment of a particle with respect to a wall changes from perpendicular to parallel on 
decreasing the particle- wall distance, see Eq. ()6.1|) . It would be interesting to check our 
predictions with simulations or real experiments. 

The simple and general forms ()3.6|) and ()3.8|) of the density of the conformal stress tensor 
in a wedge and outside a dumbbell or lens, with eigenvectors shown in Fig. 2, follow from 
combining symmetries of the two geometries, see the end of Sec. III. We show in Appendix 
C that the sum (j4.2j) of canonical tensor and improvement term is a conformal tensor, while 
the separate terms are not and have densities with a more complicated form as discussed in 
Appendix B. 



16 



Acknowledgments 



We thank T.W. Burkhardt for useful discussions. 



APPENDIX A: CONTINUITY EQUATION IN THE WEDGE 

Here we show that the f2- independence of the prefactor Tij in ()3.6|) follows from the 
continuity equation. For convenience we choose Cartesian axes perpendicular to the edge in 
the = and = 7r/2 half planes and denote them by indices v and w, respectively, so 
that 

i^v , r^) = p {cosQ , sinf2) , 

iu^:^\u^:^) = (-sinn,cosfi), (Al) 

and 

(pdQ/dfv , pdil/dfw j j — sinfi , cosfi j {A2) 
dp/df-y , dp/df-yj j \ cosfi , sinfi j 
Substituting (j3.6|) with Tjj(a) — > Tij{a,VL) into the two continuity equationsi^ 

d (^wv) wedge / 9f y + <9(T^^^,)„edge/'9'r^ = 0, (-^3) 



using ()A1|) . and calculating the derivatives by means of the chain rule and ()A2|) . one finds 
that the p-derivatives of the prefactor p~'^ cancel the f2-derivatives of the eigenvector u*^"-'. 
Only the (^-derivatives of r remain, and the left hand sides of the first and second Eq. ()A3|) 
are given by {dr / dVL){d — 1)/ p'^^^ multiplied by sinfi and — cosfi, respectively. Thus each of 
the two equations implies that r is independent of Vl. 



APPENDIX B: WEDGE WITH a = 7r/2 

The propagator (<^i2)wedge for the Gaussian model in a wedge with Dirichlet boundary 
conditions and a = tt/q, with g an integer, can be obtained by the method of images. 
Besides the half space (a = vr) we consider here the simplest case of a wedge with a = tt/2 
in which there are three images. With the notation 
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for the components of 

~ (^eij f'xij ^yi) i Z = 1, 2 (-B2) 

perpendicular to the edge along Cartesian axes in the two boundary half planes Q = —tt/4 
(i.e. Y = 0, X > 0) and Q = tt/4 (i.e. X = 0,Y > 0), the propagator at the critical point 
has the form ()4.5|1 with 

(^i2)buik = S, [el, + (Xi - X2)' + (Fi - Y,)']-(^-'y' , (B3) 

and 

^(^12)wedge = S, { -[el, + (Xi + X^y + (Fi - 

+ [el, + (Xi + X2)2 + (Fi + Y2)^]-(''-^y^ 

-[el, + (Xi - X2)2 + (Fi + } . (B4) 

Here el, = (fci — ^02)^, and 5*^ is defined in ()4.ip . 

Letting ri — > r, r2 — r on the right hand side of ()B4|) leads to the boundary-induced 
profile 

-£ip/2y-''V2 = (<l>2(r))„edge2''-V'5d = -X^-'' - Y^-" + p"-" (B5) 
of the energy density in (jHH), which due to 

X = r^ = p{cosn - sin^])/ V2 , Y = fy = /)(cos^] + sinfi) /V2 (B6) 

is in accordance with the general form (|3.1|) of a scalar density with bulk exponent Xo = d—2. 
It diverges on approaching the boundary planes Q = ±7r/4. 

For the boundary- induced densities of the canonical stress tensor ()4.3|) and the 
'improvement '-term ()4.4|) . 

((TirHr))„edge , (-4A(i-))wcdgc) = {d ' 2)5,2"'^ , -J^y) , (B7) 
Eq. (jBH) yields 

= {d- 2)r-^ + dp-'^-^x^ -id- l)p-'^ 
r,, = dp-^'-^XY 

%pep = id-2) (X-'^ + -id- 3)p-'^ (B8) 
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and 



-J.. = -{d - 2)Y-'' - [did - 2)/id - l)]p-''-^X^ + id- 2)p-^ 
-J.y = -[d{d-2)/{d-l)]p-''-'XY 

-J'e/3e/3 = - (ci - 2) (X"'^ + F"'^) + [(d - 2) V (d - l)]p-'=' . (B9) 

Here e/? with /3 = 3, ..,d runs over the d — 2 Cartesian directions of the edge-subspace. Tyy 
and —J^yy follow from the above expressions for %:x and —J^xx on exchanging X and Y. All 
other components of (T^'^'^'^^) and (— /) vanish by symmetry. 

In the sum 7^^ ~ Jk\ all contributions proportional to X~'^ and Y~'^ cancel, and the 
density of the stress tensor T in (j4.2j) has the simple form of (j3.6j) with 

r(a = 7r/2) = S""^ ^ (BIO) 

implying 

h^lL'' = -2{TT/8Y/yT{d/2) (Bll) 

if one uses (nrTH|l with V/L = l/y/2 and the forms of Bt below (lOll and 5^ in (jUT}. 
However, the densities of T^"^^^^ and / separately have nonvanishing traces, with 

d 

%x + ryy + Y.%p,p = -((v$)2)2'^-V5d 

/3=3 

= (rf - 2)[(rf - + y-'*) - (rf - 2)p-'^] , (B12) 

and a more complicated form. For example, their diagonal components parallel to the edge, 
and the trace ()B12|) . depend not only on p but, due to the terms oc + Y~'^, also on 
fi, and diverge on approaching the boundaries X = or F = 0. There is a corresponding 
divergence of the diagonal components of T*^'^'^'^) and / parallel to a Dirichlet boundary wall 
for the half space (a = vr). Moreover, u^") of Eq. and Fig. 2 is, in general, not 

an eigenvector of (T^'^^^^), due to the terms proportional to Y^'^ and X~'^ in T^x and Tyy, 
respectively. In the continuity equation satisfied by (T^'^^^^), these terms do not contribute. 

APPENDIX C: CONFORMAL STRESS DENSITY 

Here we verify that in the Gaussian model the boundary- induced profile (T), with T the 
sum fj4.2|) of the canonical stress tensor and the improvement term, is a conformal tensor. 
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For a conformal transformation f' r' which leads from a geometry G' to a geometry G' , 
we use the transformation law 

mr^'i) W2))g' = m) b{r',)Y'-'y' Smr[) $(r',))^, (CI) 

for the propagator 5 (<!>$) with bulk contribution subtracted as in ()4.5|) . The dilatation 
factor 6(r') of the transformation is defined below Eq. ()3.9p . Consider the boundary- 
induced density (Tfc/(r'))G', where T^i follows from ()4.2|) on replacing (k, A,r) by (/c,/,r'). 
As in Appendix B each of its contributions follows from the subtracted propagator on the 
left hand side of i^Cl^ by appropriate differentiations, where for the terms in —{Iki)G' and 
^j,^can)^^^ the arguments r'^ and r2 are set equal to r' before and after differentiating, re- 
spectively. Following the same steps on the right hand side of ()C1|) . one finds that in the 
sum {Tki)G' = {TkT^^)G' — {Iki)G' the contributions involving r'-derivatives of h that come 
from calculating T^'^'^°'> cancel the contributions involving r'-derivatives of dr' / dv' and h that 
come from calculating —I and one is left with the transformation formula 

W))o' = E b'-^^') (II) (^) (^^A(rO)c' (C2) 

of a conformal tensor. In general T^^^^^ and I separately do not satisfy a transformation 
law of the form of ()C2|) . It is instructive to explicitly check the above arguments for the 
inversion r' = L^r'/r'^, where 

^ = 6(r')X.,.(r') , Kr') = -, (C3) 

with X defined in Eq. ()3.12p . and where ()C2|) reduces to the transformation law of Eqs. 
(2.36) and (2.37) in Ref. U 

The transformation formula ()(^2j) for the stress tensor density is also expected to apply 
beyond the Gaussian model. Given in geometry G' an eigenvector of (T(r'))^, which points 
along a certain distance vector dr' and belongs to an eigenvalue t, Eq. ()(y2j) implies in 
geometry G' an eigenvector of {T{r'))G' which points along the corresponding distance vector 
dr' and belongs to the eigenvalue b'^(r')i. 
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APPENDIX D: STRESS AMPLITUDE FOR ARBITRARY OPENING ANGLE 

It is convenient to calculate T{a) from ()3.6p with k = A equal to an edge direction e/3. 
Using (j4.2|) - (j4.4j) and the property 



2 ((V$)')wedge = A ($')wedge 

of the vanishing trace of — J)„edge, this yields 



T{a) = p'iT^p^p) = 



((*e,)') 



wedge 



A{d- 1 



-A($2) 



wedge 



where 



(Dl) 



(D2) 



(D3) 



Since the boundary induced profile 



wedge 



lim 



ri— >r,r2— >r 



[($e/3(ri)$e/3(r2))wedge - ($e/3(ri)$e/3(r2))bulk] (D4) 



is, apart from a factor Sd, analytio^i in d, the calculation proceeds similar to ()4.11|l - ()4.17j] : 
We consider d < 0, where 



($e/3(fl)$e;3(r2))bulk = (d - 2)[\h " f 2 1' " ^(^le/S " hcfir]]^ " hr''' (D5) 



(D6) 



does not contribute on the right hand side of ()D4jl . and find for = 

1 



P 



with 



where 



(($eMi-))^)wedgeb=0 = '^^^""'^^"^ 



(rf/2))^- 



(D7) 



Edit) = t('^-2)/2(i (D8) 

Here we have rewritten (($c/3)^)wcdge as J2 pii^eisy) wedge /{d — 2) and used Eqs. ()4.7|1 and 
(1121) with 2A = d. The integral (|D7jl converges for -2n/a < d < 0. 
Writing as the sum of 



Jo 



1 



1 - t2-/c 



and 



JP = dtEd{t)r/"k{t) 
Jo 



(D9) 



(DIO) 
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with 



k{t) 



a 1 



1 



1 - ^ 1 + — ( — ) [2(1 - t) + (1 - tf\ 



27r 1 - t ' 2 V" iTi) 24 V a 27r 
the continuation to = 3 leads to 



/27r 



27r\^ «■ 
— -1+6— 
V a / 27r. 



yielding 



1 



(($e/3(r))')wcdgc|n=op'' = [Ji^)(a,3) + Ji')(a,3)] 



in c? = 3 where J^^[ol^ 3) follows from ()D9j) on replacing by H3 = t^/^(l — t) 
In the second contribution to r in (ID2I). 



p'^A($2)^edge = -V25,{92^+(rf-2)2^} 



can, for i7 = and (i = 3, be expressed by Cq and 62 in ()4.17|) . ()4.18p . so that 

r(a) = -^[J(i)(a,3) + Jf («,3)] + ^_eo(a)[e2(a) + 1] . 
«7r 16v27r 
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a 





7r/2 


TT 


27r 


a 




2-i/2in2 


2-1/2 _ 2-2 


2-3/2 


2-V2,r-i 






= 0.490 


= 0.457 


= 0.354 


= 0.225 






(0.569) 


(0.508) 










7C(3)/(64V2) 


(4 + 2-V2)/(64V2) 




-5/(967r\/2) 






= 0.0930 


= 0.0520 





= -0.0117 






(0.0561) 


(0.0304) 










-7rC(3)/4 


-7r2-V2 




1/6 






= -0.944 


= -0.555 





= 0.166 






(-0.630) 


(-0.346) 









TABLE I: Amplitudes for dumbbells with a = 0, 7r/2, the sphere (a = tt), and the disk (a = 
27r) in three spatial dimensions. Numbers in brackets denote amplitudes for prolate ellipsoids 
circumscribing the dumbbells. The values correspond to the crosses and circles in Figs. 3-5. 
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FIG. 1: Conformal mapping of a wedge onto a dumbbell or lens. 



26 



/ 



\ 



/ 



/ 



..V 

ir 

■■/., 



T-/.. 

I 

.I...V 



FIG. 2: The normalized eigenvectors u(")(r2) (left) and u(^)(r) (right) of the stress tensor density 
in a wedge and outside a dumbbell or lens. The vector field u(^) is determined by Eq. P-lOf) . The 
symbols for lines and points correspond to those in Fig. 1. 
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FIG. 3: The leading amplitude ai in H2.4() for particles with Dirichlet boundary conditions in the 
Gaussian model. Crosses denote results for dumbbells (0 < a/vr < 1) and lenses (1 < a/ir < 2). 
Circles denote results for circumscribing prolate and oblate ellipsoids, see the paragraph containing 
Eqs. 1)4. 19(1 and (|4.2U|) . The oblate ellipsoid and the lens coincide for a/vr = 1 and 2, where they 
become a sphere and a disk, respectively. 
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FIG. 4: The anisotropy amplitude h\ in (|2.4|1 . 
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FIG. 5: The anisotropy amplitude 62 in (|^.4|) . 
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